Relativistic particle in a three-dimensional box 
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We generalize the work of Alberto, Fiolhais and Gil and solve the problem of a Dirac particle 
confined in a 3-dimensional box. The non-relativistic and ultra-relativistic limits are considered 
and it is shown that the size of the box determines how relativistic the low-lying states are. The 
consequences for the density of states of a relativistic fermion gas are briefly discussed. 
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The problem of a particle confined in an one-dimensional infinite square well potential lies at the heart of non- 
relativistic quantum mechanics, being the simplest problem that illustrates how the wave nature of bound particles 
I— 1 1 implies that their energy is quantized. The generalization of this problem to three dimensions is used for the statistical 
i-C . description of a fermion gas and is the starting point for some many-body theories of fermions like the Thomas-Fermi 
model of the atom. 

The relativistic formulation and solution of the problem of a spin- 1/2 fermion with mass m confined in a one- 
dimensional square well potential was done by Alberto, Fiolhais and Gil [1]. In this paper we review their approach 
! and generalize it to a 3-dimensional square well potential. The problem of a relativistic particle confined in infinite 
square well potential is traditionally not dealt with in the textbooks of Relativistic Quantum Mechanics, even in the 
'— ' ^ most comprehensive ones such as the one by Greiner [2] . 

_ -, ■ Recently, in the context of quantum gravity phenomenology, applications of this result were presented. In particular, 
^ ' by applying the Generalized Uncertainty Principle (GUP) to a particle confined in a three-dimensional box, it is shown 
^vq that the length of the box must be quantized in terms of a fundamental length, e.g. Planck length, and this indicates 
^\ that the nature of space may be fundamentally grainy [3] . 

T-H , In the present analysis, we choose an ansatz for the spinor inside the box which is consistent with the non-relativistic 
CO ■ problem and derive a transcendental equation for the wave numbers allowed. We obtain the non-relativistic and ultra- 
psj ' relativistic limits of this equation and its solutions, showing also how they are related to the ratio between the size 
C3 ' of the containing box and the Compton wavelength h/{mc) of the fermion. Finally, we discuss briefly the differences 
regarding the density of states between the usual non-relativistic fermion gas and the relativistic fermion gas. 
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II. SOLUTION OF DIRAC EQUATION WITH A ONE DIMENSIONAL INFINITE SQUARE WELL 



, The Dirac equation is written as ^ 



Hip = {cd ■ p + jSmc^^ ilj{r) 

= Ei^ir) (1) 



where (« = 1, 2, 3) and /3 are the Dirac matrices, for which we use the following representation 

' P \ -I 



°"'\/3=(i°) (2) 
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^ In this section, wc follow the formulation of [1] . 
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where cTi arc the PauU matrices and / is the 2-dimcnsional unit matrix. It is evident that in 1-spatial dimension and 
in the position representation, say z, the Dirac equation is given by 



The positive energy solutions read 



-ihca^-^ + /3mc^ ) ip{z) = Etp{z) . (3) 




where ra is the mass of the Dirac particle, k is the wavcnumber that satisfies the usual dispersion relation E'^ = 
{hkcY + {mc^y, r = J^^^i and x^X = 1- Note that r runs from (non-relativistic) to 1 (ultra-relativistic) , k could 
be positive (right moving) or negative (left moving), while A/' is a suitable normalization constant. 
As noted in [1], to confine a relativistic particle in a box of length L in a consistent way avoiding the Klein paradox 
(in which an increasing number of negative energy particles are excited), one may take its mass to be z-dependent as 
was done in the MIT bag model of quark confinement [4, 5] 

m{z) = M, ^ < (Region I) (4) 
= m,^<z<L (Region II) (5) 
= M, z>L (Region III) , (6) 

where m and M are constants, and taking eventually the limit M oo. This guarantees that the confinement process 
avoids the excitation of negative energy Dirac sea particles, since the plane wave energy spectrum of positive and 
negative energy solutions is separated by at least twice their mass times (? . Thus, the general form of the wavefunction 
for a bounded Dirac particle in a one dimensional box, vanishing when \z\ — )• oo, can be written (in all three regions) 
as 

( X \ 

(7) 

(8) 
(9) 




where K = sj E"^ — [Mc^Y /{he) (an imaginary wavenumber when M > E/c^) and R = hKc/{E + Mc^). Thus, in 
the limit M — > oo, ii' — > +ioo, the terms associated with A and D go to zero. Now, boundary conditions akin to those 
for the Schrodinger equation, namely ijjn = at z = and z — L will require iljn to vanish identically. Thus, they 
are disallowed. This is related to the fact that the usual boundary conditions for non-relativistic quantum mechanics 
cannot always be applied to relativistic problems, as was shown by V. Alonso et al. [7]. Instead, we require the 
outward component of the Dirac current to be zero at the boundaries (the MIT bag model). This ensures that the 
particle is indeed confined within the box [5]. 
The conserved current corresponding to Eq.(3) can be shown to be 

Jz = ^Yi' , (10) 

where 7^ = 7^ is the 2:-component of the four- vector set of 4 x 4 gamma matrices 7** = (7", 7). These matrices are 
related to the ai and /3 matrices introduced earlier by 7*^ = /3 and 7* = jSai . The vanishing of the outward component 
of the Dirac current = ipjf^ijj at a boundary is obtained by requiring that the condition ■ nip = ip holds there, 
where n is the unit four- vector normal to the space-time boundary [5], such that n^n'^ = — 1. For a static boundary 
= (0, i), where x is the outward normal to the boundary surface. Applying this to the wavefunction ipn at 2; = 
and z = L gives [1] 

and - i^azV'z/l^^/, = V'//|^=i (12) 
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respectively. Using the expression for ^j/ from (8), we get from (11) and (12) respectively 

B + C 



B-C 

which in turn yield 



= ir (13) 
= -ir , (14) 



(zr-l) = g(ir + l) (15) 

(ir + 1) = §(zr-l)e-^2'=^ . (16) 
B 



By eliminating C/B between (15) and (16) one gets 

ir — 1 



ir + 1 

so that we finally arrive to the transcendental equation 



2r hk 

tan(A;L) = ^5 = . (18) 

r"" — 1 mc 

The discrete solutions of this equation for the values of the wavenumber k give the quantized energy levels for a 
relativistic particle in the 1-dimcnsional box. Note that the solution with A: = is excluded because, from (15), 
C = —B and the upper component would be proportional to sin(A:x) and therefore be identically zero inside the box. 
It is worth mentioning at this point that Eq. (15) implies that \B\ — \C\ and this condition could have as well be 
obtained by just requiring that J^, given by (10), vanishes at the boundary, as remarked by Menon and Belyi [6]. 
Note that this condition does not imply that B = ±(7, since at least one of B or C is not a real number in this case, 
as can also be seen from Eq. (15) and the fact that r 7^ 0. Also one may note that in order to obtain the eigenvalue 
equation (18), one needs to know the ratio of the complex coefficients C and B, C/B, which we got from the MIT 
boundary condition, and not only the ratio of their moduli. 

To conclude this section, wc comment briefly on the negative energy solutions of the Dirac equation (3) in a one- 
dimensional box. These solutions can be found from the negative energy solutions of the free Dirac equation travelling 
in opposite directions in the z axis. The new boundary condition is obtained from the corresponding condition for 
positive energy by applying the charge conjugation operator. The resulting equation for the wavenumber k is identical 
to (18) and therefore its discrete solutions k„ are the same. This means that one has an overall symmetric energy 
spectrum, in which for every discrete positive level with energy E = En = ^JW'k'^ + m^c** there is one with E = —En- 
This could be expected, since in our case one has a Dirac Hamiltonian with a confining Lorentz scalar potential which 
anti-commutes with the charge conjugation operator, as does the free Dirac Hamiltonian. Since the Hamiltonian with 
the present boundary conditions is Hermitian (sc!c [6]), all wavefunctions belonging to distinct eigenvalues must be 
orthogonal to each other. Thus the positive and negative energy solutions do not mix. Similar conclusions hold for a 
two and three dimensional confining box as well. 

III. SOLUTION OF DIRAC EQUATION WITH A THREE-DIMENSIONAL INFINITE SQUARE WELL 

The generalization of the problem of the previous section involves solving the 3-dimensional Dirac equation with a 
position-dependent mass of the form 

{m re V 
M r^V 

where V is defined by the set of points with coordinates (xi, X2, xz) such that Q <x\ < ii, < xs < L2, < xs < Z/3. 
As before we take the limit M — >■ cx). The potential (19) can be written in a more compact way as 

3 

m{xi,X'2, X3)=m+{M - m) ]J[6'(a;i - Li) + e{-Xi)] (20) 
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where 9{x) is the step function. For the solutions inside V we expect to have combinations of free Dirac spinors with 
positive energy, which have the general form 




(21) 



where k is the wave vector and r = h\k\c/{E + mc^). 

In the non-relativistic case (i.e., for the Schrodinger equation with free particles confined in the volume V), since 
the potential (20) is separable, one can use the product ansatz for the wave function 



(22) 



and thus turn the three-dimensional problem into a set of three independent one-dimensional problems, the total 
energy being just the sum of the one-dimensional energies. In the relativistic case, however, one cannot use such 
an ansatz because of the spinor structure of the wavefunctions. Nevertheless, we may use that feature of the non- 
relativistic wavefunction as a guide to find the right combination of spinors (21). Indeed, we may require that in the 
non-relativistic limit, when the lower component of the spinor vanishes, the space part of the remaining 2-component 
spinor be of the form (22). A wavefunction that meets this requirement is 



TL=i n-=i(S,e*'=^-^ +Q(-l)^ 



(23) 



where an overall normalization has been set to unity. It can be easily shown that the above is a superposition, with 
appropriate products of the coefficients Bi and Cj, of the following 8 eigenfunctions, for all possible combinations of 
Ci [i = 1,2,3), with a = ±1 



\^ _ g» eifeiK; 



X 



(24) 



Each one of these eigenfunctions is of type (21), and they represent plane waves travelling in the 8 directions 
(±fci, ±fc2, ifcs) all with the same momentum magnitude p ~ hk ~ h^Jk\ + fc| + /cg. 

Again, we impose the MIT bag boundary conditions ±iPaitp = ip ,1 = 1,2,3, with the -I- and — signs corresponding 
to a;( = and Xi = Li respectively, ensuring vanishing flux through all six boundaries. First, we write the above 
boundary condition for any Xi, for the wavefunction given in Eq.(23). This yields the following 2-component equation 



m=l 



n^=i (Bje^*^^^ +Cj(-l)*i-e-»'=^^^)mi kmCJr, 



X 



Kj^j _|_ (J.^ 



(25) 



By equating the upper components of the spinors on each side of this equation and invoke the arbitrariness of x one 
gets the matrix relation 



H [ n [Bje''"'''' + Cj{-lY^"^e-'''i''i) rk, 

m=l j=l 



(26) 



In the same fashion, equating the lower components of the spinors on each side of Eq.(25) and right multiply the 
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resulting matrix equation by ±i ai one gets 

3 3 

m=l j=l 

3 

= Jl (Sje*'=^^^- + Qe-^'=^^^) 7 . (27) 



Adding these two equations one gets 

3 



3 3 

m—l j — 1 
3 

= ±2i[ JJ [Bje'^^^'i +Cj{-lf''e-'''i''^) rfc;]/ (28) 

where the relation {c;, (j„,} = 2di„,I was used. Noticing that for j ^ I the terms in the products arc the same in both 
sides of the equation, one can divide both sides by Yi^^i {BjC^^^^^ + CjC"**^^^^) and obtain, for each I, 

Bie''"'" + C;e-''='^' = ±i (Sic''^'^' - ^6"''='^') rk . (29) 

Note that the boundary condition (25) gave rise to 3 independent condition pairs, one for each spatial dimension. 
Eq.(29) yields, at x; = and Xi = Li respectively 

Bi + Ci=i{Bi-Ci)rki (30) 

and 

Bie''"^' + Cie-'^'^' = -i {Bie'''^^' - C;e-''='^') rk . (31) 

Comparing Eqs.(30) and (31) with Eqs.(13) and (14), we see that we can write, in a similar way as in the previous 
section, 

'I^LZI = e'hL, (32) 
irki + 1 

and thus 

2rki _ 2{E + mc^)hki 

' ' - r^kf - 1 ~ h^{kf-k^)c-2mc{E + mc^) ' ^^^^ 

These arc a set of three coupled transcendental equations for fci , fc2 , ^3 that need to be solved to find the energy 
eigenvalues. Again, the solutions with A;; = are excluded. In the non-relativistic limit, E ~ mc^, hki/{mc) = e; <C 1 
so that 

kiLi ^ arctan(— e;) ^ n;7r (34) 

where ri; = 1, ... so that we recover the well-known quantization conditions for a non-relativistic fermion gas confined 
in a cubic box. 

On the other hand, in the ultra-relativistic limit, when E ~ hkc » mc^ and hki/{mc) ^ 1, one gets 

kiLi ~ arctan ( —a^— \ ~ {m — l/2)7r (35) 
\ mc J 

where a is a coefficient whose value depends on relative magnitude between ki and the other components (i.e., a = 1 

if fc( >fcr, V ^l). 
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The size of the box plays a fundamental role in determining the relativistic behavior of the solutions. Indeed, one 

may write 

hL, = ^ (36) 
mc Lc 

where Lc = h/{mc) is the Compton wavelength. Given that the righthand side of Eq.(33) is always negative (we 
take every ki to be positive) and therefore its first solution such that 7r/2 < kiLi < tt, we see that when Li » Lq and 
Li <IC Lc we get the non-relativistic and ultra-relativistic limits, respectively. To illustrate this point, we write the 
eigenvalue equation in terms of ratio Lc/Li. When ki is much bigger than the other components, Eq.(33) reads 

Lc 

tana;; = —xi (37) 
Li 

where x; = The graphical solution of this equation is shown in Fig. 1. 




FIG. 1: Graphical solution of equation (37) for three values of the ratio Li/Lc- The dotted vertical corresponds to xi = ■k/2. 

From Fig. 1, one can check that indeed the size of box determines the relativistic nature of the solutions and also 
that the first non-relativistic solution is ki ~ tt/L; and the first ultra-relativistic solution is ki ~ 'nl[2Li). Of course, 
the degree of relativity increases for the higher energy solutions. Note that if the condition of ki being much bigger 
than the other components is relaxed, this would amount to a small change in the coefficient of Xi in the right-end 
side of Eq.(37), provided that ki = maxj(fcj), so that the previous conclusions would still hold. Note that all these 
conclusions assume that one has always a non-zero mass for the fermions. 

In Fig. 2, it is depicted the energy spectra corresponding to the solutions of Eqs. (33) for L; = Lc/IQ, Li = Lc, 
Li = lOLc, considering a cubic box L\ = L2 = L^- For convenience, the values plotted are of the logarithm of the 
scaled kinetic energy E/{m,c?) — 1. The levels plotted correspond to the first 27 levels for each value of L;. However, 
due to the symmetry of equations (33) with respect to the interchange of the ki 's among them, some of the levels are 
degenerate (besides, of course, the spin degeneracy). If fci ^ ^2 7^ ^3, we have a 3! = 6 -fold degeneracy, while when 
ki 7^ fc;/ = kill , {I, I', I"} = permutations of {1,2,3}, one has a 3-fold degeneracy. This reduces the 27 levels to 10 
levels of distinct energies. 

This degeneracy is completely analogous to the one found in the non-relativistic case, in which the kinetic energy 
is given hy E — mc^ = h^ir^ / {2mL'^){nl +nl + rig) , ni, n2, ns = 1, . . ., L being the size of the cubic box. In fact, one 
can use the non-relativistic quantum numbers ni to classify the quantum states in the present case. This is presented 
in Table I, which contains the quantum numbers, the level degeneracy, the values of kiLi/ir and corresponding energy 
values for the first 6 distinct energy states. 

One can check from the values presented in the table that indeed, for small values of Li/Lc, h approaches (n; — 
1/2)'k/Li while for higher values they approach niiT/Li, the non-relativistic value, as suggested by Fig. 1. 



FIG. 2: Scaled kinetic energy spectrum of the first 27 solutions of equation (33) for three values of the ratio Li/Lc- 



{ni,n'2,nz) 


dcgcn. 


(fei,fe2,A;3)(xL(/7r) 


Li/Lc 


E/{m(?) 






(0.674129, 0.674129, 0.674129) 


0.1 


36.6957 


(1,1,1) 


1 


(0.730735, 0.730735, 0.730735) 


1 


4.10004 






(0.914156, 0.914156, 0.914156) 


10 


1.11689 






(0.761157, 0.761157, 1.5664) 


0.1 


59.718 


(1,1,2) 


3 


(0.789821, 0.789821, 1.61153) 


1 


6.24063 






(0.917935, 0.917935, 1.8383) 


10 


1.22469 






(0.800534, 1.62894, 1.62894) 


0.1 


76.6236 


(1,2,2) 


3 


(0.820262, 1.66176, 1.66176) 


1 


7.88349 






(0.921162, 1.84449, 1.84449) 


10 


1.32488 






(0.819801, 0.819801, 2.53383) 


0.1 


87.5453 


(1,1,3) 


3 


(0.835499, 0.835499, 2.56592) 


1 


8.93086 






(0.923098, 0.923098, 2.77709) 


10 


1.38902 






(1.66969, 1.66969, 1.66969) 


0.1 


90.8601 


(2,2,2) 


1 


(1.69565, 1.69565, 1.69565) 


1 


9.28075 






(1.84989, 1.84989, 1.84989) 


10 


1.41889 






(0.838015, 1.69214, 2.57123) 


0.1 


100.225 


(1,2,3) 


6 


(0.850724, 1.71438, 2.59869) 


1 


10.1883 






(0.92567, 1.85317, 2.7841) 


10 


1.47937 



TABLE I: Non-relativistic quantum numbers, degeneracies, values of h in units of tt/L; and scaled energies for the first 6 
distinct energy solutions of equation (33). For each level are presented the values for Li/Lc = 0.1, 1, 10. 

Finally, we will make some considerations regarding density of states of relativistic fermion gases. From Fig. 1 and 
Table I one sees that as Li gets bigger each allowed value for ki is separated by ir/Li as in the non-relativistic fermion 
gas. For fermion gases confined in large volnmcs with many particles, as considered for instance in Statistical Physics, 
the Dirac gas in a 3-D box behaves as a non-relativistic gas. For sizes V ~ L^, for which the fermions are relativistic, 
the density, measured by the separation of ki values in units of n/Li, tends to increase, as can be seen from Fig. 1 and 
Table I, since each solution is separated by less than tt/Li. For smaller volumes, we checked that we recover again 
the non-relativistic density of states (NRDS) when the quantum numbers changing arc relatively high. In general, we 
can say that, for higher energy levels, the density of states tends to be the same as the NRDS. On the other hand, 
for small quantum numbers, in relativistic (and ultra-relativistic) conditions, we get higher densities than the NRDS. 
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IV. CONCLUSIONS 



In this paper, we have solved the problem of a relativistic spin- 1/2 particle in a three-dimensional square box using 

the Dirac equation. Wc studied both the non-rclativistic and ultra-rclativistic limits and showed that these are related 
to the size of the box. The scale for gauging the relativity of the solution is the Compton wavelength, such that free 
fermions confined in boxes with sizes many times the Compton wavelength behave as non-relativistic particles. We 
solved the three coupled transcendental equations which give the energy eigenvalues of the Dirac particle in a 3D box 
and used the non-rclativistic quantum numbers of quantum particle in a box to classify the corresponding states. We 
also showed that the density of states is in general higher than the non-relativistic one, but for higher energy states it 
will tend to be equal to non-relativistic density of states, in which the allowed values for the wavenumber components 
ki, I = 1, 2, 3 are separated by tt/L;, where L; is the respective box dimension. 
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